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Abstract. In this work we study the theory of extended quasidilaton massive gravity together 
with the presence of matter helds. After discussing the homogeneous and isotropic fully 
dynamical background equations, which governs the exact expansion history of the universe, 
we consider small cosmological perturbations around these general FLRW solutions. The 
stability of tensor, vector and scalar perturbations on top of these general background solutions 
give rise to slightly different constraints on the parameters of the theory than those obtained 
in the approximative assumption of the late-time asymptotic form of the expansion history, 
which does not correspond to our current epoch. This opens up the possibility of stable FLRW 
solutions to be compared with current data on cosmic expansion with the restricted parameter 
space based on theoretical ground. 


Contents 


1 Introduction 


1 


2 Extended quasidilaton dRGT massive gravity 


3 


3 Background evolution 


5 


4 Stability of the perturbations 


6 

7 

9 

10 


4.1 Tensor perturbations 

4.2 Vector perturbations 

4.3 Scalar perturbations 


5 Conclusions 


13 


1 Introduction 

It has been a long time since we know that the universe is expanding, which means that if 
we look at two distant galaxies, the separation between them increases in time. The rate 
at which the galaxies are getting apart is determined by the content of the universe and the 
current consensus is that ordinary matter (baryons, electrons ... etc) only represents about 
hve percent of the content of the universe. A larger fraction, about 27%, corresponds to 
what is known as dark matter. This is a component that behaves like a dust fluid but does 
not interact (or at least very weakly) with photons. Its existence is however inferred by its 
gravitational effects, ranging from galactic scales to the universe horizon. It is also crucial 
for the formation of the structures that we observe today. Finally, maybe the most intriguing 
component is what is usually called dark energy, which amounts to about 70% of the universe 
content and it is responsible for the accelerated expansion of the universe. This accelerated 
expansion of the universe was hrst inferred by measurements of type la supernovae. These 
supernovae are standarizable candles, that can be used to measure cosmological distances. 
One of the pillars of the standard model of cosmology is the Cosmological Principle, which 
states that the universe is homogeneous and isotropic on large scales. These symmetries im¬ 
pose that the metric of the universe must be of the Friedmann-Lemaitre-Robertson-Walker 
(FLRW) form, where the scale factor a{t) determines physical distances. Thus, by measuring 
the distance to supernovae at different redshifts, we can have a measure of the expansion rate 
of the universe. In General Relativity, from Einstein equations one can obtain the acceleration 
equation (second derivative of the scale factor). It is immediately clear that in order to have 
an accelerated expansion, the universe must be dominated by a fluid satisfying the condition 
a > 0, or in other words the pressure of the fluid must fulhl the condition P < —\p- Thus, 
this condition is not satished by the usual fluids like dust or radiation. 
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The simplest way to achieve acceleration is the introduction of a cosmological constant 
A. This gives a contribution in the form of a perfect fluid with Pa = —pa, and, therefore, pro¬ 
vides acceleration. This has become the concordance model or standard model of cosmology, 
called ACDM. It is able to £t most cosmological observations with a high accuracy. We also 
have measurements of large scale structures where we measure the distribution of galaxies in 
the universe, and ACDM also gives good hts to current data. This type of measuremets will 
be the next generation of cosmological observations. Although the cosmological constant is in 
good agreement with current observations, it has some problems from a theoretical viewpoint. 
First of all, if A is a true constant of nature belonging to the gravitational sector, then the 
gravitational action has two dimensionful constants which differ by many orders of magnitude 
and this makes it unnatural. On the other hand, the cosmological constant will also receive 
radiative corrections from quantum helds, and those corrections are again many orders of 
magnitude larger than the observed value, so we need to hnely tune the bare value to be in 
agreement with observations. The unnatural problems of the cosmological constant and its 
instability against quantum corrections motivated searching for alternatives to a cosmological 
constant in order to explain the cosmic acceleration. Such alternatives can be broadly divided 
into two categories. On one hand, models that rely on the introduction of scalar helds, vector 
helds, new exotic huids, etc as the cause of the accelerated expansion. On the other hand, 
there are also models that resort to modifying GR at large scales, arguing that the acceler¬ 
ated expansion might be due to a break down of Einstein’s General Relativity on cosmological 
scales. 

One interesting way of modifying gravity in the infra-red (IR) is massive gravity. The 
simplest linear massive gravity theory is the pioneering model proposed by Fierz and Pauli 
[1]. The nonlinear covariant completion of the Fierz-Pauli theory with the correct degrees 
of freedom was formulated by de Rham, Gabadadze and Tolley (dRGT) [2, 3]. The poten¬ 
tial is tuned in a way that guarantees the absence of the Boulware-Deser (BD) ghost [4]. 
This relative tuning is radiatively stable in the decoupling limit of the theory due to the an¬ 
tisymmetric structure of the potential interactions [5] (Galileon interactions share this nice 
property as well [5-9]). The radiative stability does not survive beyond the decoupling limit, 
however the introduced BD ghost is harmless since its mass is beyond the cut-oh scale of the 
theory [10]. Additionally, the quantum corrections arising from matter loops might dictate 
the form of allowed classical interactions with matter helds [11-13]. The phenomenology of 
the dRGT theory is reach. Its potential impact on cosmology is a very interesting question 
[14-34]. The construction of a realistic and stable cosmology has witnessed many promising 
attempts, even though most of them were doomed to a failure very quickly. The hrst failure 
was the construction of hat FLRW solution [15]. Similarly, even if it was possible to hud 
self-accelerating open FLRW solutions, these were suhering from a nonlinear ghost instability 
[16]. Attempts to promote the reference metric to de Sitter or FLRW did not help since they 
face the unavoidable appearance of Higuchi ghost [25]. Possible ways out of these difhculties 
consist on breaking the FLRW symmetries [15, 18, 19, 22, 24, 35-37] or adding additional 
new degrees of freedom [38-40]. One interesting way of extending the original massive gravity 
theory is the inclusion of the kinetic term for the reference metric [41]. Another promising 
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way is the coupling matter fields through an effective composite metric [11, 42, 43]. 

In this work we will go along the line of the presence of additional degrees of freedom. 
To be precise, we will focus on the extension ’quasidilaton’ in form of an additional scalar 
field with a specific coupling to the massive graviton [39]. The Lagrangian is constructed in 
such a way that it is invariant under the quasidilaton global symmetry. Even if the original 
version of the quasidilaton fails to provide stable self-accelerating solutions, there is a promis¬ 
ing extension of it proposed by de Felice and Mukohyama [40] through the addition of a new 
coupling constant. Interestingly enough, the background dynamics is independent of this new 
coupling, however for the stability of the perturbations it does play a crucial role. In the stan¬ 
dard formulation of massive gravity, phenomenological consistency reqnires a small graviton 
mass, of the same order as the Hubble expansion rate today m ^ Hq. As a consequence the 
effective mass of tensor perturbations would be also constrained to be of the same order. In 
the quasidilaton scenario the effective mass of gravitational waves has an additional depen¬ 
dence on the quasidilaton parameter and hence it can be made much larger than the Hubble 
parameter. This might provide an explanation for the large-angle suppression of power in the 
microwave background to be explored [44]. 

The scalar perturbations of the extended quasidilaton with matter fields were studied in 
[45]. In their analysis, they assumed the late-time asymptotic form of the expansion history 
to evaluate the stability of scalar perturbations in the presence of matter fields. However, the 
late-time asymptotic solntion does not correctly describe our current epoch, which is just in the 
transition between matter domination and accelerating expansion. In this work, we consider 
the full dynamical equations of motion and study the tensor, vector and scalar perturbations 
on top of these general background equations. By doing so, we obtain new stability constraints 
which highly depend on the evolution of background quantities, which were set to zero or to 
some constant values in the asymptotic solution. We first review the extended quasidilaton 
generalization of dRGT theory and setup our Lagrangian and conventions in Section 2. After 
calculating the fully dynamical backgronnd eqnations of motion on FLRW space-time in Sec¬ 
tion 3, we pay special attention to the stability of the tensor, vector and scalar perturbations 
in the presence of the matter fields in Section 4, where we obtain new highly non-trivial stabil¬ 
ity constraints on the parameters of the theory. Finally, we snmmarise our results in Section 5. 

Thronghout the paper, we will work with the metric signature convention (—, -|-, -h, +) 
and define the reduced Planck mass as Mpi = 1/ VSttG. Furthermore, the traces of some 
rank-2 tensors in Section 2 are denoted by [...], for example /C^„ = [/C], = [/C^] = 

(/C^,)2, = [/C3] = (/C^,)3 etc. 

2 Extended quasidilaton dRGT massive gravity 

In this section we will first review the interactions in the theory of extended quasidilaton 
massive gravity and setup the framework in which we will perform our analysis of cosmological 
perturbations. Our starting point is the action for extended quasidilaton massive gravity and 
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the matter action where the ordinary matter helds still couple minimally to the physical metric 
g and an additional scalar held a, the quasidilaton held, comes in the potential interactions 
in a very specihc way such that the absence of the BD ghost is maintained [40] 

S = J d^x [ \/—g (^R[9] ~ 2A — ^2 dij,ad^a + 2m?{hl2 + 03^3 + + -^matter] , 

( 2 . 1 ) 

where the potential interactions are given by [2, 3] 

«2[K] = 1 ([K]^ - ri). 

MsM = i([K]=‘-3[/C]|K^] + 2[/C=>]), 

" 3 IK] = i (|K1* - 6 [K]=[/C"] + 3[/Cy + 8|/C1[K’] - 6 |/C‘]) , (2.2) 

The tensor K, is the building block of dRGT theory [3] with a highly non-trivial structure in 
form of a square root. In the extended quasidilaton massive gravity this tensor is promoted 
to ^ 

= (2.3) 

The presence of the potential interactions in 2.1 breaks the diheomorphism invariance com¬ 
pletely. Naturally, one can restore it introducing four Stiickelberg helds 5“ which promote the 
Minkowski reference metric to the space-time tensor 

f,. ^ Vabd.S^d^S^ - . (2.4) 

As you can see the standard form of the reference metric —?■ riabd^S°‘dyS’^ is promoted to 

2.4. Note that this disformal transformation to the hducial metric introduces a new parameter. 
The building block tensor /C is invariant under the global symmetry a ^ a + (Tq and iS“ —?■ 
g-cro/Mpi^a parameter a^j in the disformal factor was introduced in order to render 

the self-accelerating late-time asymptotic solutions stable [40]. For sake of simplicity we will 
consider a matter held 0 that only couples to the dynamical metric g with a generic kinetic 
term of the form 

-^matter = V~9 P{^4>) ) (2-5) 

where = —g^''d^(j)du4> and the energy density pm, the pressure Pm and the sound speed 
Cm of the matter held can be written as 

PM^2P(X^,)'X^,-P(X^,), PM^i’(W), 4 ^ 2P(xI''x!i PiX.y • 

where prime corresponds to derivative with respect to the argument. In our setup in the 
Lagrangian 2.1, we neglected the tadpole contribution lAi to the mass term by setting ai = 0, 
while absorbing the Wq term into the bare cosmological constant A. 
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3 Background evolution 


The original formulation of massive gravity does not accommodate a simple realisation of 
viable cosmology due to a no-go result for flat FLRW solutions [15]. This arises from the 
constraint equation imposed by the Stiickelberg field equation of motion, which demands the 
scale factor to be constant. In order to avoid this no-go result, extensions to the original for¬ 
mulation were considered in the literature. Promising attempts are quasidilaton and extended 
quasidilaton formulation of the dRGT theory where an additional scalar field a was added. 
We would like to study perturbations on top of FLRW backgrounds of the latter case in the 
presence of matter helds coupled to the dynamical metric. Thus, we will start with the Ansatz 
for the dynamical metric as the homogeneous and isotropic flat FLRW metric 


ds^g = —N'^dt^ -I- o?5ijdx^dx^ , (3.1) 

while the non-dynamical metric is the pull-back of the Minkowski metric in the Stiickelberg 
field space to the physical space-time, which we parametrize as 


ds^f = f^ydx^dx^ 


n^dt^ -I- dijdx^dx^ 





dt^ + dijdx^dx^. 


(3.2) 


where the lapse of the / metric is = p + ^ order to be compatible with 

the homogeneous and isotropic Ansatz for the metric, we also assume that our matter held 
only depends on time 0 = 0(t). Furthermore, we introduce the following quantities for our 
convenience 


(3.3) 

(3.4) 

A)q;3 + (1 — A)^q;4, (3-5) 

where A denotes A = and H is the expansion rate of the physical g metric, while r 

encodes the speed of light propagating in the / metric in the units of the one propagating in 
the g metric. The action for the extended quasidilaton (2.1) in the mini-superspace becomes 


H = 


r = 


a 

na 


N ’ 

J = 3 + 3(1 


F 

V 


P, 


= I dta'^N^ - A - 3H'^ -Aa + m^Jr{A -l)A+^ + 


ua 




where we made use of the following dehnitions 

Aa = ™^(A - 1) |J + (4 - l)(a3(Al - 1) - 3)] 


, (3.6) 


(3,7) 


We are now in possession of all the important quantities to be to compute the background 
equations of motion. We do this by varying the action (3.6) with respect to N, a, (p, a and /. 
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Of course as usual, the resulting system of equations of motion contains a redundant equation 
giving by the contracted Bianchi identity, 


d SS a6S f 6S & SS cj) SS 

rnsN ivfe nJJ nI^ nI^ 


(3.8) 


First of all, we start with the Friedmann equation, which can be obtained by varying the 
action (3.6) with respect to the lapse N, giving rise to 



(3.9) 


As next we vary the mini-superspace action (3.6) with respect to the scale factor a and combine 
the resulting equation with the Friedmann equation, which amounts to 


2H (1 - r)kA pM + Pm 


(3.10) 


N 3HN - 3&/Mpi 


The equation of motion for the matter held is just the standard conservation equation 


!^ + iH(pM + PM) = 0. 


(3.11) 


Last but not least, the variation with respect to the Stiickelberg held results in 



(3.12) 


The remaining equation of motion for the a held reads 

m^MpiN^A{3{r - l)A{-2 + a^A - 1)) + J(-3 + r(-l + 4A))) 

= uj (3HN^& + Na- N&) (3.13) 


As you can see, the dynamical background equations highly depend on the quantities J, A and 


similarly on a and the dynamics of the matter held (j). In [40], a late-time attractor solution 
for the system in absence of matter helds was reported. In this late-time asymptotic solution 
the background quantities are either zero or constants. To be precise, J = 0 and H, A and 
r are constants. With these assumptions, our background equations of course coincide with 
the ones in [40]. 

4 Stability of the perturbations 

The central goal of this work is to study the stability conditions for the perturbations around 
the full dynamical background equations from the previous section. In diherence to [45], we 
will assume J ^ 0, A ^ 0 and a ^ Mp\NH, i.e. we will not assume the late-time asymptotic 
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form of the expansion history, since we are just at the transition between matter domination 
and accelerating expansion making the late-time asymptotic solution not suitable to describe 
our current epoch. 

Let us begin with the following perturbations for the dynamical metric gfj,iy 


6goo = -2 , 

6goi = Na {diB + Bi) , 
2 


5gij = 2 Siji> + (^didj - y <9^4^ E + d^iEj) + hi 


(4.1) 


where it is understood that all the metric perturbations depend on time and space. Note also 
that 6^^hij = d'^hij = d^Ei = d^Bi = 0. Furthermore, we do not consider any perturbations 
for the Stiickelberg helds, which hxes the gauge freedom completely. We perturb the matter 
field (j) and the quasidilaton a as well 


0 = (j){t) -h MpiScj), 

a = a{t) -|- Mpi6a . (4.2) 

Naively counted, we encounter twelve degrees of freedom (dof) in the action 2.1. Two of these 
dos are traceless symmetric spatial tensor helds (hij). In addition, we have four divergence- 
free spatial vector helds {B^, Ei). The remaining six dof come in form of scalars ($, B, ip, E, 
Scj), a). In the following we will study the tensor, vector and scalar perturbations separately 
and establish stability conditions above the dynamical background equations. 


4.1 Tensor perturbations 

Before starting the computation of the tensor perturbations, let us hrst decompose the tensor 
held hij in Fourier modes with respect to the spatial coordinates. We can perform this 
decomposition since our background metric is homogeneous with no spatial curvature 

d^k 

(2^)3/2 ^P,fc(^) • (4-3) 



We hrst insert our Ansatz 4.1 for the metric perturbations into the Lagrangian 2.1 with the 
helds decomposed in Fourier modes and then we use the fully dynamical background equations 
from the previous section. By doing that, the action quadratic in the tensor perturbations 
results in 


where the mass of the tensor perturbations is given by the abbreviation 





( 2 ) 




tensor 


dPk dt N oP 


k 


rrin 


k 



p-—^(J(-l -FrA) -F A(-(2 -F a3){-2 + r) - {1 + 2 ao + r)A + aovA^)). 


(4.5) 
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The tensor perturbations on top of general background solutions (4.4) have already the right 
sign for the kinetic term. Similarly, they do not exhibit gradient instabilities either. Further¬ 
more, in order to avoid tachyonic instabilities we have to impose m\ > 0. This condition is 
fulhlled for the following three cases 


0 < r < A~^, 
r > A~^, 
r = A~^, 


J > J {J < J), 0 < A < 1 (A > 1) 

J < J {J > J), 0 < A < 1 (A > 1) 

j 4(1 -|- tt3) < 2 -|- tt3, > 0 


(4.6) 


where J = " 3 )) _ mass term can be also expressed in terms of 

the dynamics of the matter held (T..etc by replacing through its expression given by the 
equation of motion of the dilaton a held 3.13. 




m? (3(r — 1)‘^A^ — JA{—3 + r + 2rA){r{2A — 1) — 1)) 
3(kl-l)(r-l) 
uj{rA + r — 2)>V 
^3Mp,Af2(r - 1)(/1 - 1) 


(4,7) 


where we introduced the shortcut notation >V = SHN^a — N& + Na. On the self-accelerating 
(SA) late-time asymptotic background with J = 0, A = 0, d = NMpiH, H = 0 and with the 
gauge choice = 1, the above mass term coincides with equation (38) of de Felice et al in 


[40] 


m 


2 _ 

T,SA — 


m'^irsA - 1)^A|^ + uHlj^{rsA{AsA + 1) - 2) 
(^^ SA — l)(As^ — 1) 


(4.8) 


where tsa is obtained from the equation of motion for the dilaton 


rsA = l+ u:HlJ{m^AsA{oi-3,{AsA - 1) - 2)) (4.9) 


and Hsa from the Friedmann equation = Api^^(3 — a;/2)“k On this background Hsa-, 
Asa and tsa are constants. The stability of the self-accelerating background solutions requires 
0 < cj < 6 , and tsa > 0 for Asa > 1 or 1 < tsa < '^sa for 0 < Asa < 1 where [44] 


rsA 


2 

1 + 2\.sa 


OJ 

Asa{^ + Asa) 



I) + - 1)^ 


(4.10) 


Furthermore, for simplicity we have assumed that the sound speed of the tensor perturbations 
given in 4.4 is unity, i.e. = 1. Therefore the dispersion relation is given by^ 


f = ^^rnl (4.11) 

where f is the frequency of tensor mode oscillations. Inflation generates super-horizon gravi¬ 
tational waves (with -C H^) that stay constant outside the horizon with an absolute value 


^The dispersion relation in pure GR, fcR = k'^/a^, simply obtains an additional mass term. 














given by the amplitude of the mode at initial time (i.e. \hk\ = A[k) with A[k) = H*/ (Mpi/c^/^) 
where H* is the Hubble expansion rate at horizon exit in a typical slow-roll inflation) . As 
the universe evolves (i.e. decreases), these modes at some point re-enters the horizon 

(with ~ iJ^) at time a^e and dually end up as sub-horizon modes (with ^ H'^). In 
pure GR, the horizon crossing happens at = /c^^/(\/2fc) for momenta larger than keq 
and = (/c^^/(\/2/c))^ for smaller momenta respectively, where a^q is the time of matter- 
radiation equality. After the horizon re-entrance the gravitational waves start oscillating with 
frequency f as a WKB oscillator. In pure GR, the power spectrum of sub-horizon 

gravitational waves today {t = to) is given by with k = Oofo- In massive 

gravity with the additional mass term the power spectrum becomes instead [46] 


V 


fo 


fo - "if,0 ^ 


2 l^fcl 


(4.12) 


with /c = ooY fo — Q. 
versus GR is given by^ 


Thus, the relative change of the power spectrum in massive gravity 


_ f fk^k'Ppr^k 


(4.13) 


with the primordial power spectrum Vpr^k = ‘ik^A^/. In the late-time asymptotic regime, 
we expect that the mass term evolves into the constant value mx —t mx^sA given by equation 
(4.8). However, we need the full solution of the background equations in order to see how mx 
evolves in time and to be able to compute its value today and also during matter domination. 
It would be very interesting to solve 4.13 numerically and study the impact of the additional 
dilation held (for the case in the pure standard massive gravity see [46]). In the standard pure 
massive gravity theory the mass of the graviton is constrained to be of the same order as the 
Hubble expansion rate today m ^ Hq. Therefore this severely restricts mx to be of the same 
order mx ^ Hq. In the quasi dilation case this can be avoided. Even if m ~ oiie can have 
mx S> Hq due to the presence of the additional degree of freedom. A natural outcome from 
this might be the large-angle suppression of power in the microwave background [44]. This 
together with numerical solutions and comparison to current data will be studied elsewhere. 


4.2 Vector perturbations 

We would like now pay our attention to the stability conditions of the vector perturbations. 
As we did for the tensor perturbation, we hrst decompose the vector modes Ei and Bi in 
Fourier modes 


/ d^k 


B, = 


ri^k 

(4.14) 


As next we expand the Lagrangian 2.1 to second order in the vector perturbations. Note that 
not all of the vector helds are dynamical, indeed the vector helds Bi do not have any kinetic 

^ For a detail discussion see [46]. 


9 










terms. We can therefore compute the equation of motion with respect to Bi and integrate 
them out 


o _ k'^a{l + r) _ • __ 

2N{k'^ + k?r+ 2m?a^A{J — [2 + a‘i)A +a'iA^)) 

After using this expression, the quadratic action in the vector perturbations becomes 


(4.15) 




( 2 ) 


vector 


16 


d^k dt k'^a^N 




miE*rEl 

t,k k 


where we introduced for convenience the shortcut notation 

2m?a?A{J + A{—2 + q; 3 (A — 1))) 


my = 


N‘^{k‘^ + k'^r + 2m?a^A{J — (2 + 03 )^ + 03 ^^)) 


(4.16) 


(4.17) 


For the stability of vector perturbations on top of a generic background solution we have to 
impose the right sign for the kinetic and gradient terms. For the absence of ghost instability, 
we require that the kinetic term has the right sign. This is the case if we impose rriy > 0. 
This on the other hand requires r > 0, A > 0 and J > A{2 + 0 ( 3(1 — A)) or fc^(l + r) + 
2a^m?A{J + A(—2 + (A — 1 ) 0 ( 3 )) < 0. Similarly as for the tensor perturbations the absence of 
gradient and tachyonic instability requires rri^ > 0, which are the same constraints as in 4.6. 
Again we can replace the dependence on 0(3 through the other parameters of the theory using 
the equation of motion for the dilaton held. In this case the prefactor in front of the kinetic 
term becomes 


my = 


1 

iW 


3k‘^MpiN{r^ - 1 ) 


My\N^{-3k‘^ + r(3Fr - 3m‘^a?J{A - 1)A)) + 2uja?W 


(4,18) 


Recall that W = SHN^a — Na + Nd. On the self-accelerating background with J 
0, d = NMp\H, H = 0 and with the gauge choice N = 1, the expression for my 
simply 

_2 2^0^ Hi 

= m{2ooami^ + k^{rl^-l)) 
which coincides with the expression of de Felice et al [40]. 


= 0,A = 
becomes 

(4.19) 


4.3 Scalar perturbations 

Now we will be concentrating on the stability of the scalar perturbations in the extended 
quasi-dilaton massive gravity model with matter held. As we mentioned above, six degrees of 
freedom appear in form of scalar helds {ip, 6a, 6(j), E, B, ^ ). We hrst expand the action 2.1 
to quadratic order in the scalar perturbations and introduce their Fourier modes 


$ = 
Ip = 
6cP= , 


dP^k 

(27r)3/2 

d^k 


ik-x 


iPj^ (t)A"'" 


(27r)3/2 
d^k 


(27r)3/2 


+ c.c. 


+ c.c. 


+ c.c. 


B = 


E = 


6a = 


d^k 

(27r)3/2 

d^k 

(27r)3/2 

" d^k 

(27r)3/2 


B^ {t)E’^-^ + c.c. 
% + c.c. 


6aj: {t)E^'^ + c.c (4.20) 
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The corresponding kinetic matrix (and the Hessian matrix) contains two vanishing eigenvalnes, 
signalling the existence of two constraint eqnations which make two ont of the six scalar helds 
not propagating 


/C 




/-6 0 0 0 0 0\ 

0 u! 0 0 0 0 

0 0 2P{X^y + 4P{X^)"X^ 0 0 0 

00 0 k^/6 0 0 

00 0 000 

\ 0 0 0 0 0 0 / 


(4.21) 


Since the qnadratic action does not have any kinetic term for the scalar helds <h and B, we 
can compnte their eqnations of motion in order to obtain the corresponding two constraint 
equations 

+ eM^HNii - eMytiSi - + 6P(X4,)'iS4'( ( 4 . 22 ) 


and similarly 

+ 3{2eMl,aB^HN + 2eMl,N'iP^ + ^2(_ 3^Pi^Q^ag _ 

(H — 1) 

1 SA/f^ N^AOih- 

+QM^,HiPj^ + —i + N\-Mp,ua6aj: + - 2Mpi0P(X^)'% 

+24>^<i>^P{X^y) + 40=^(-Mpi%^ + 4>,^0)P(X^)"))) = 0 (4.23) 


We can solve them for B and $. For clarity, we have dehned the quantity 

Q = + H) + (Ayi + m‘^{A — 1)^)A. 


(4.24) 


After plugging back the solutions for B and <F, the resulting action depends only on the 
remaining four scalar helds ("0, 5a, 6(j), E). The kinetic matrix of these four remaining scalar 
helds still has a vanishing determinant, meaning that there is still one more constraint that 
can be used to integrate out one of the scalar helds. In fact, the sub-kinematic matrix for the 
helds ("0, 5a, 50) has already a vanishing determinant and it has one zero eigenvalue Ai = 0 
and two non-vanishing eigenvalues A 2 and A 3 

/O 0 0\ 

= 0 A 2 0 (4.25) 

yo 0 As/ 

The corresponding eigenvectors to the three eigenvalues Ai, A 2 and A 3 are respectively given 
as 

/MpiHN/^\ 

Ul = ( (t/0 1 , 
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We omit the cumbersome expressions for the eigenvalues and eigenvectors here. In the UV 
limit for instance we manage to write down the eigenvalues in a manageable way 


where 

I = + 2Ml^H^N^P{X^y + N^oja"^ + 2N^^'^P{X^)' 

+AMl,H^N^^‘^P{X^)" + 40^P(X^)" 
q = + ^'^){N‘^P{X^)' + 2^‘^P{X^)'') 


(4.27) 


(4.28) 


Before moving on, at this stage it is useful to remind the reader that these perturbations are 
performed around the general background equations including the matter helds. Comparing 
our analysis with [45], it is clear that the difference comes fro the fact that the authors in [45] 
did not perform their analysis using the full dynamical background equations. The dynamical 
background evolution dictates that a = Mpi {^N + 3 ^ and not a = Mp\HN. Dynamical 

terms of the form A/A are missing in [45]. At this stage, it is not clear to us how to justify 
to use the late-time asymptotic form of the expansion history to evaluate the perturbations 
for our current epoch. The full dynamical background equations would need to be solved 
numerically and the perturbations should be performed about these exact solutions. 


The presence of the zero eigenvalue in 4.25 indicates that we can hnd a constraint equa¬ 
tion to eliminate one of the degrees of freedom, the would-be Boulware-Deser ghost. After 
integrating it out, the resulting action can be expressed in the following form 

5'IcLr = ^ j d^kdta^ kti + tA Un-lA Uti-lA MYi^ , (4.29) 

where If denotes the three physical propagating degrees of freedom If = 

The expressions for the kinetic, gradient matrices ..etc are very cumbersome around general 
background solutions after integrating out the would-be Boulware-Deser ghost, even in their 
UV and IR limits, which makes it very hard to write them down in a handy way. However, 
the for us important conditions here are the ones coming from the ghost absence. For this, 
we have to guarantee that the kinetic matrix is positive dehnite. This on the other hand 
imposes conditions on the determinant of the matrix det iF > 0 and parts of it det JF 33 > 0 , 
-^ 22-^33 — > 0. By requiring these conditions on the kinetic matrix in the UV and IR 

limit together amounts to^ 

0 < m < 6 and A^ < ^ (4.30) 

\ A I 

^In the self accelerating late-time asymptotic solution the condition is the same as reported in [40]. 
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As you can see, this condition differs from the one derived in [45]. Since they were assuming 
the late-time asymptotic form of the expansion history, terms proportional to J, A and & were 
missing in their analysis. However, one needs to work with the full solution of the dynamical 
background equations to describe the evolution of the universe. The reported instability in 
[45] was probably just due to the missing terms A/A based on the late-time approximated 
solution. It would be very interesting to study numerical solutions of the general background 
equations and check the consequences of the stability conditions on the parameters. The exact 
numerical solution and comparison to the current data on cosmic expansion is being studied 
elsewhere. 

5 Conclusions 

In this work, we have considered the extended quasidilaton massive gravity theory, which 
includes an additional scalar degree of freedom on top of the massive graviton. This extended 
version of the quasidilaton has an additional parameter, which was introduced to give rise to 
stable self-accelerating solution at late times. On top of these self-accelerating background it 
was shown that the perturbations of the tensor, vector and scalar perturbations were stable 
after imposing the necessary conditions on the parameters of the theory [40]. Here we extend 
their stability analysis to general background including the presence of matter fields. 

After discussing the general dynamical background eqnations, we studied the stability 
conditions on the tensor, vector and scalar perturbations. The tensor perturbations have 
automatically the right sign for the kinetic and gradient term without imposing further con¬ 
straints. The only constraint coming from the tensor perturbations is the positivity of the 
mass term > 0. From this condition we were able to put constraints on the background 
quantities J, r, ...etc and on the free parameters of the theory. From the vector perturbations 
we obtain an additional constraint coming from the coefficient in front of the kinetic term, 
m\r > 0 in order to avoid ghost instability. This condition puts further constraint on J and 
the parameter together with the condition r > 0. The absence of gradient instability of 
the vector perturbations does not give any further constraint since it is the same constraint 
as for the absence of tachyonic instability of tensor perturbations. Last but not least, the 
stability of the scalar pertnrbations not only imposes constraints on the parameter uj and 
but also on the dynamical quantities A. Again for the stability of the scalar perturbations, 
the parameter plays a crucial role. An important step in future works would be to study 
the dynamical system of cosmological solutions by using phase map analysis and explore the 
critical points of the cosmological equations as well as numerical solutions. 
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